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Within a one particle approximation of the Dirac equation 
we investigate a defect system in a quantum wire. We demon- 
strate that by minimally coupling a laser field of frequency u> 
to such an impurity system, one may generate harmonics of 
multiples of the incoming frequency. In a double defect sys- 
tem one may use the distance between the defects in order to 
tune the cut-off frequency. 

PACS numbers: 32.80 Rm, 73.40.-c 



I. INTRODUCTION 

Since the early nineties, the generation of high-order 
harmonics of a strong, low-frequency laser field has at- 
tracted a lot of attention in the atomic physics commu- 
nity, see e.g. [1,2] for a fairly recent review. Indeed, high- 
harmonic generation has opened a wide range of possibil- 
ities for obtaining high-frequency coherent sources, con- 
verting infrared input radiation of frequency u> into light 
in the extreme ultraviolet regime whose frequencies are 
multiples of u> (even up to order ~ 1000, see e.g. [2] for 
a recent review). 

In gases, composed of atoms or small molecules, this 
phenomenon is well- understood and, to some extent, even 
controllable in the sense that the frequency of the highest 
harmonic, the so-called "cut-off", can be tuned as well 
as the intensities of particular groups of harmonics. In 
more complex systems, however, as for instance solids, 
or larger molecules, high-harmonic generation is still an 
open problem. This is due to the fact that, until a few 
years ago, such systems were expected not to survive the 
strong laser fields involved. However, nowadays, with the 
advent of ultrashort pulses, there exist solid-state materi- 
als whose damage threshold is beyond the required inten- 
sities of 10 14 W/cm 2 [3]. As a direct consequence, there 
is an increasing interest in such materials as potential 
sources for high- harmonics. In fact, several groups are 
currently investigating this phenomenon in systems such 
as thin crystals [4,5], carbon nanotubes [6], or organic 
molecules [7,8]. 

Prototype solid-state devices are quantum wires, which 
nowadays do not only serve as theoretical test ground, 
but may even be studied experimentally, e.g. [9]. With 
regard to the previously outlined problematic it is of 
great interest to investigate how such devices interact 
with laser light. In particular, the question of whether 



such systems are suitable high-harmonic sources has not 
been dealt with up to now. This is the central question 
which we shall be answering positively in this manuscript. 

A useful particularity of quantum wires is that they 
are the physical realization of models involving only one 
spatial dimension. Considering theories in one space di- 
mension has the further virtue that it allowed for the 
development of various powerful non-perturbative tech- 
niques, which exploit the integrability of the models. For 
instance, one may compute the conductance of quantum 
wires in this fashion, see e.g. [10] and references therein. 
This also suggests that an exact relativistic treatment of 
the electrons in the wire in the presence of the laser field 
is possible. Therefore such a treatment constitutes an 
advantage over most approaches adopted in the litera- 
ture, in the context of atoms in strong laser fields, which 
involve a series of approximations. The first, and most 
general approximation performed in this context is that 
a proper field-theoretical treatment is usually not taken 
into account. Furthermore, when dealing with a rela- 
tivistic situation, the Dirac equation is solved mostly by 
numerical methods, which by themselves involve a series 
of approximation and exhibit a good convergence only 
for high-frequency fields [11]. Further approximations 
include the expansion [12] of the Dirac equation for the 
weakly relativistic case, the solution of the Klein-Gordon 
equation [13,14] or of the relativistic Schrodinger equa- 
tion, or purely classical treatments [15]. For some recent 
general reviews see e.g. [16]. 

In this paper, we solve the Dirac equation for an elec- 
tron in a quantum wire subject to an external laser field, 
including a single and multiple defects. We investigate 
in various different regimes the radiation emitted by an 
electron in this system in connection with the transmis- 
sion at the defects and demonstrate that the generation 
of harmonics is possible in such systems. 



II. DEFECT SYSTEMS IN LASER FIELDS 

Since the work of Weyl [17], one knows that matter 
may be coupled to light by means of a local gauge trans- 
formation, which reflects itself in the usual minimal cou- 



pling prescription, i.e. 



ieA^, with 



being 



the vector gauge potential. For a free Fermion with mass 
m this yields to the Lagrangian density 



C A = tpin^df, -m+ ej^A^ip . 



(1) 
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As common conventions 1 , we abbreviated here ip = 1/^7° 
and use the following realization for the gamma matrices 



7° = 
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(2) 



In the absence of the laser field the equations of motion 
for the free Fermion may be solved by the well-known 
Fourier decomposition 



il>Ux,t) 



de 



(3) 



We parameterize the momentum as common through the 
rapidity 9 by p° = mj cosh#, pj = rrij sinh# and denote 
the anti-particle (positron) of the Fermion (electron) j 
by J . For the Weyl spinors we employ the normalization 



Uj (0) = -i^ Vj (9) 




(4) 



and the creation and annihilation operators a,(#), a\{6) 
for a particle with rapidity 9 obey the usual 
fermionic anti-commutation relations {aj(#i), ^(^2)} = 
0, {a;(0i),at(0 2 )} = 2tt%(5(0i - 9 2 ). When the laser 
field is switched on, we can solve the equation of motion 
associated to (1) 



(i^d^ - to + ej^A^tp = 
by a Gordon- Volkov type solution [18] 







ij)f{x, t) = cxp 



ie / dsAi(s,t) 



exp 



ie J dsA (x,s) 



i> f j(x,t). 



(5) 

(6) 
(7) 



Using now a linearly polarized laser field along the di- 
rection of the wire, the vector potential can typically be 
taken in the dipole approximation to be a superposition 
of monochromatic light with frequency u), i.e., 



A(t) :=A 1 (t) = - f dsAo(s) 
x Jo 



(8) 



dsf(s) cos(ws) (9) 



with f(t) being an arbitrary enveloping function equal 
to zero for t < and t > t, such that r denotes the 
pulse length. In the following we will always take f(t) = 



1 We adopt here relativistic units l = c = ft = m~ e 2 137 as 
mostly used in the particle physics context rather than atomic 
units 1 — e — h — mKi c/137 useful in atomic physics. 



0(i)O(r — t) , with 0(x) being the Heavyside unit step 
function. The second equality in (8), A (x,t) = xA(t), 
follows from the fact that we have to solve (6) and (7). 

One comment is due with regard to the validity of the 
dipole approximation in this context. It consists usually 
in neglecting the spatial dependence of the laser field, 
which is justified when xiv < c = 1, where x is a repre- 
sentative scale of the problem considered. In the context 
of atomic physics this is typically the Bohr radius. In 
the problem investigated here, this approximation has to 
hold over the full spatial range in which the Fermion fol- 
lows the electric field. We can estimate this classically, in 
which case the maximal amplitude is eE /u 2 and there- 
fore the following constraint has to hold 



UJ 



4U P < 1 



(10) 



for the dipole approximation to be valid. Due to the fact 
that x is a function of u>, we have now a lower bound on 
the frequency rather than an upper one as is more com- 
mon in the context of atomic physics. We have also intro- 
duced here the ponderomotive energy U p for monochro- 
matic light, that is the average kinetic energy transferred 
from the laser field to the electron. 

The solutions to the equations of motion of the free 
system and the one which includes the laser field are then 
related by a factor similar to the gauge transformation 
from the length to the velocity gauge 



ipf{x, t) = cxp [ixeA(t)\ ipj(x) . 



(11) 



In an analogous fashion one may use the same minimal 
coupling procedure also to couple in addition the laser 
field to the defect. One has to invoke the equation of 
motion in order to carry this out, since as in [10], we 
assume here also that the defect is linear in the fields 
-ip and tp. The Lagrangian density for a complex free 
Fermion ip with t defects T> a (ip, ip, A^) of type a at the 
position x n subjected to a laser field then reads 



Cad = £a + J2 ^> A ^ ■ 



(12) 



71=1 



Considering for simplicity first the case of a single defect 
situated at x = 0, the solution to the equation of motion 
resulting from (12) is taken to be of the form 

ipf{x,t) = Q(x)il>£ + (x,t) + e{-x)il>f ! _(x,t). (13) 

This means we distinguish here by notation the solutions 
(11) on the left and right of the defect, tpf_(x,t) and 
^f+{x,t), respectively. This is also reflected in the cor- 
responding creation and annihilation operators a,j y +{9), 
dj-{9), etc. One may then proceed according to stan- 
dard potential scattering theory and notes that these 
functions are not independent of each other. Substitution 
of (13) into the equation of motion yields the constraints 
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i7 1 (^+(M)-V£- (*,*))!*= 



x=0 

(14) 



These restrictions (14) serve to determine the transmis- 
sion and reflection amplitudes. The extension to multiple 
defects, that is having equations of the type (14) for each 
defect situated at position x — x n , follows in an analo- 
gous straightforward manner. 



1 r 

T(n,e,E ,w,T) = - dt\T(e,E Q ,uj,t)\ 2 cos(nt), (19) 

T Jo 

1 r 

n(n,6,E ,u,T) = - dt\R(6,E ,uj,t)\ 2 cos{nt). (20) 
T Jo 

When parity is preserved for the reflection amplitudes, 
we have \T\ 2 + \R\ 2 = 1, and it suffices to consider T in 
the following. 



1. Type I defects 



A. Transmission and Reflection amplitudes 

Substituting the Fourier decomposition (3), together 
with the free Fermion solution minimally coupled to a 
laser field (11) into the constraint (14), one can determine 
the reflection and transmission amplitudes from the left 
to the right, R and T, respectively. They are defined in 
an obvious manner as 

a jt -{0) = R^ajA-O) +T j (0)a j , + (0). (15) 

When parity invariance is broken, the corresponding am- 
plitudes from the right to the left do not have to be iden- 
tical and are defined as 

H+(-0) = f 3 (9)H-(-9) + Rj(0)a jt+ m ■ ( 16 ) 

In this way the laser field parameters Eq and u> will be 
quantities on which the i?'s and T's depend upon at a 
particular time t. When iterating these equations one 
obtains the corresponding expressions for multiple defect 
systems, e.g. placing for instance the defect V a left from 
one obtains the well known expressions 



Tf{9) = 



l- R?(6)Rf(0)' 



Rf(6) = R?(6) 



R-(0)T?(0)f*(9) 
l-Rf(6)R?(6) 



(17) 
(18) 



From our previous comment on the validity of the dipolc 
approximation it is clear that the distance between the 
two defects introduces a new scale in the system, which 
has to be constraint as y < In addition, to justify 

that the multiple defect system can be treated effectively 
as a single one, requires that the sum of the y's is much 
smaller than the length of the wire. Similar expressions, 
which we will not need in what follows, hold for the parity 
reversed situation and for more defects, see e.g. [10] and 
references therein. 

In addition with regard to the application of high har- 
monic generation, we shall be interested in the spectrum 
of frequencies which are filtered out by the defect while 
the laser pulse is non-zero. The Fourier transforms of the 
reflection and transmission probabilities provide exactly 
this information 



Taking the laser field in form of monochromatic light in 
the dipole approximation (9), we may naturally assume 
that the transmission probability for some particular de- 
fects can be expanded as 

oo 

\T I (6,U p ,oj,t)\ 2 = J2 t 2k(0)^U p ) k S m 2k (ujt). (21) 

k=0 

We shall refer to defects which admit such an expansion 
as "type I defects". Assuming that the coefficients t2k(9) 
become at most 1, we have to restrict our attention to the 
regime AU P < 1 in order for this expansion to be meaning- 
ful for all t. Note that this is no further limitation, since 
it is precisely the same constraint as already encountered 
for the validity of the dipole approximation (10). The 
first coefficient is always the transmission probability for 
vanishing laser field, that is t (0) = \T(9, E = 0)| 2 . The 
functional dependence of (21) will turn out to hold for 
various explicit defects considered below. Based on this 
equation, we compute for such type of defect 



T/(fi,(9,C^,w,r) = 2 



(2k)\(U p ) k sm(rtyt 2 k(0) 

to mnf=iP 2 -(fi/2w) 2 ] ' 



(22) 



It is clear from this expression that type I defects will 
preferably let even multiples of the basic frequency u> 
pass, whose amplitudes will depend on the coefficients 
*2fc(0)- When we choose the pulse length to be integer 
cycles, i.e. r = 27rn/w = nT for n £ Z, the expression in 
(22) reduces even further. The values at even multiples 
of the basic frequency are simply 



T!{2nuj,e,U p ) = {-l) r 



£M*) (u P ) k 

k=0 



which becomes independent of the pulse length r. Notice 
also that the dependence on E and u occurs in the com- 
bination of the ponderomotive energy U p . Further state- 
ments require the precise form of the coefficients t2k{0) 
and can only be made with regard to a more concrete 
form of the defect. 



2. Type II defects 

Clearly, not all defects are of the form (21) and we have 
to consider also expansions of the type 
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\T n (6,E /e,u;,t)\ 2 = E ^W^2^cosPM)sin 2fe M). 

k,p=0 

(24) 

We shall refer to defects which admit such an expansion 
as "type II defects" . In this case we obtain 



oo p 



T n (n,e,E /e,u;,T) = J2 J2 
k,p=0 1=0 

( 



p 



ft sin(rf2) 
(-l) i+1 ru; 2 + 2fe 



(25) 



(2k + 2iy.t 2 2 l(0) (2k + 2l)k? k +1 (e)E 



k-\-l 

. n i(2<z) 2 - (S) 2 ] n [(2?-i) 2 -s 2 : 

\ 9 =0 ~ 9 =1 



E, 



2k+2p 




We observe from this expression that type II defects will 
filter out all multiples of u>. For the pulse being once 
again of integer cycle length, this reduces to 



k,p=0 1=0 



(26) 



2k + 21 
k + 1 — n 



and 



T H ((2n-l)u;AEo/e)= £ ]T(-l) i+ " )2/ 

k,p=0 1=0 

, k+ pfp\ (2k + 2l)\(2n-l)E 2 Q p+1 



x (U P T 



I J {l + k-n + l)\(l + n + k)\ 



(27) 



which are again independent of r. We observe that in 
this case we can not combine the Eq and u> into a U p . 
The analytical expressions presented in this section will 
not only serve as a benchmark for our analytical com- 
putation, but can be used in addition to extract various 
structural information as we sec below. 



B. One particle approximation 

In spite of the fact that we are dealing with a quantum 
field theory, it is known that a one particle approxima- 
tion to the Dirac equation is very useful and physically 
sensible when the external forces vary only slowly on a 
scale of a few Compton wavelengths, see e.g. [19]. We 
may therefore define the spinor wavefunctions 



*j>,o(a:,t) : = il>f{x,t) 



e -ipj-x 



Vf(*,t)t 



27T 2 n" 



aUO) 



27rp° 



i(0) (28) 



v 3 (6)K (29) 



With the help of these functions we obtain then for the 
defect system 



4-W 



V 2w Pi 

e(-x) [*i, u , 9 (a;,t) + *i, u ,- e (x,t)R*(6)] + 
e(x)T*(9) \^i^e{x,t) + ^ itU ^ {x,t)R*{-e) 



(30) 



and the same function with u — * v. Since this function 
resembles a free wave it can not be normalized properly 
and we have to localize the wave in form of a wave packet 
by multiplying with an additional function, g(p,t) in (3) 
and its counterpart g(x,t) in (30), typically a Gaufiian. 
Then for the function <f>f u e {x, t) = g(x, t)^f u6 (x, t), we 
can achieve that II* II = 1. 



III. HARMONIC GENERATION 

As mentioned above, harmonic generation has at- 
tracted large attention in the atomic physics community 
in recent years. It is the non-linear response of a dipolc 
moment, in general atomic, to an external laser field. 
Here we want to investigate whether such responses also 
exist for defect systems. We carry out our treatment rel- 
ativistically. The time dependent dipolc moment is given 

by 

Xj , u ,e(t) = (9f, Ui g(x, tfx*f iUte (x, t)) (31) 

such that the emission spectrum is the absolute value of 
the Fourier transform of the dipole moment 



f 

Jo 



dtXj tU fi{t) expifit 



(32) 



We localize now the wave packet in a region much smaller 
than the classical estimate for the maximal amplitude 
the electron will acquire when following the laser field. 
We achieve this with a Gaufiian g(x,t) = cxp(— x 2 /A), 
where A <C eE /u; 2 . Placing the defect at the origin and 
neglecting its extension, the computation of (31) with 
(30) then boils down to the evaluation of 



Xi, u ,e(t) 



1 + \R46)\ 2 + 2Re- 



2ix sinh 6 



cosh 6 



Ri(0) 



2lie i[x+eA{tW 



dx + 



oo 

/ 



x^O)? l + |i?,(0)| 2 + 2Re-' 



cosh 9 



dx. (33) 



2lTe ± [x +eA(t)}i 



The expressions for the i?'s and T's depend of course 
on the form of the defect and further generic statements 
can not be made at this point. We therefore turn to a 
concrete example. 
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IV. THE ENERGY OPERATOR DEFECT 

Localizing sharply the energy operator e(x) = g'tptp(x), 
with g being a coupling constant, yields a defect which 
has been studied extensively in the absence of a laser 
field. It should be noted that this is only a particular 
example and one may also consider other type of defects 
in a analogous fashion, see [10] and references therein. 
One of the virtues of this defect is that it is real, thus 
preserving parity invariance. 

Coupling the vector potential minimally to this type 
of defect yields 

V A D(i>,i>,A li )=g^(l + e/m'fA li )^, (34) 

by invoking the equation of motion. 

A. Transmission and Reflection amplitudes 

We are now interested in determining the reflection and 
transmission amplitudes associated to this defect by the 
potential scattering method as outlined in section II A. 
Taking from now on m = 1, we compute for the various 
reflection amplitudes 

Ri(0, g, A/e, y) = Ri(9, g, -A/e, -y) = 

R T (9, g, A/e, -y) = R T (9, g, -A/e, y) = (35) 

[yA- cosh e}e~ 2iysinh 6 
[1 - yAcosh9] - i{ + 1 + A 2 - y 2 A 2 } sinh (9 ' 

We denoted the differentiation with respect to time by a 
dot. The transmission amplitudes turn out to be 



Ti(0, g, A/e, y) = f x (9, g, -A/e, -y) = 
T T (6, g, -A/e, y) = f T (0, g, A/e, -y) = 



(36) 



1 - y 2 A 2 + (A - |) 2 



sinh# 



|[1 - yA cosh 9} - i[4f + 1 + A 2 - y 2 A 2 } sinh (9 

Locating the defect at y = 0, the derivative of A does not 
appear anymore explicitly in (35) and (36), such that it is 
clear that this defect is of type I and admits an expansion 
of the form (21). With the explicit expressions (35) and 
(36) at hand, we can determine all the coefficients t 2 k{9) 
in (21) analytically. For this purpose let us first bring the 
transmission amplitude into the more symmetric form 



m(9,g,A/e)\< 



& (9,g) + a 2 (9,g)A 2 +a 4 (9,g)A 4 
a (9,g) + a 2 (9,g)A 2 +a 4 (9,g)A 4: 



with 



a (9,g) = 16. 9 2 + (4 + .g 2 ) 2 sinh 2 0, 
ao(0,g) = (g 2 -4) 2 smh 2 9, 
a 2 (9,g) = 2g 2 (4 + g 2 ) S mh 2 9, 
0,4(9, g) = g 4 sinh 2 9. 



(37) 

(38) 
(39) 
(40) 
(41) 



We can now expand \T(9,g,A)\ in powers of the field 
A(t) and identify the coefficients t 2k (9,g) in (21) there- 
after. To achieve this we simply have to carry out the 
series expansion of the denominator in (37). The latter 
admits the following compact form 



= J2c 2k (9,g)A 2k , 



k=0 



a (9,g) + a 2 (9,g)A 2 + ai (9,g)A 4 

with c {9,g) = l/a (9,g) and 

C2k-2(9,g)a 2 (9,g) + c 2k - 4 (9,g)a 4 (9,g) 



c 2 k{9,g) = -- 



ao{9,g) 



for k > 0. We understand here that all coefficients c 2k 
with k < are vanishing, such that from this formula all 
the coefficients c 2k may be computed recursively. Hence, 
by comparing with the series expansion (21), we find the 
following closed formula for the coefficients t 2k {9,g) 

t 2k (9,g) = [a (9,g)-a Q (9,g)}c 2k (9,g) k > 0. (42) 

The first coefficients then simply read 

ao(9,g) 



to{9,g) 
t2(0,g) 

U{9,g) = 



ao{9,g) 
a2{9,g) 



ao(9,g) 
8 5 4 (4 



= \T(9,E = 0)\ 2 , 
[l-t {9,g)] 



9 



2 ) sinh 2 29 



(16g 2 + (4 + g 2 ) 2 sinh 2 9f 
ai(9,g) a 2 (9,g) 



a2(9, g) a (9,g) 



t2(9,g), 



(43) 
(44) 
(45) 
(46) 



and so on. It is now clear how to obtain also the higher 
terms analytically, but since they are rather cumbersome 
we do not report them here. 



B. Harmonic generation 

With the coefficients t 2k , we can compute the series 
(22) and (23) in principle to any desired order. For some 
concrete values of the laser and defect parameters the 
results of our evaluation are depicted in figure 1. 

The main observation is that the defect acts as a fil- 
ter selecting higher harmonics of even order of the laser 
frequency. Furthermore, from the zoom of the peak re- 
gions, we see that there are satellite peaks appearing near 
the main harmonics. They reduce their intensity when 
t is increased, such that with longer pulse length the 
harmonics become more and more pronounced. 

We also investigated that for different frequencies 
lu the general structure will not change. Increas- 
ing the field amplitude E , simply lifts up the whole 
plot without altering very much its overall structure. 
We support these findings in two alternative ways, ei- 
ther by computing directly (19) numerically or, more 
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instructively, by evaluating the sums (22) and (23). 




o 5 io 15 Oca 20 25 30 35 40 

Figure 1: Absolute value squared of the Fourier transform of 
the transmission probability for a single defect with Eo = 2.0, 
g = 3.5, 9 = 1.2, u = 0.2. 

Let us now carry out a similar analysis for a double 
defect system. We place one of the defects at x = and 
the other at x = y. The distance emerges now as a new 
scale in the system and note from our comment on the 
validity of the dipole approximation that is restricted as 
y < uj~ 1 . In addition y has to be much smaller than the 
total length of the wire. The double defect amplitudes are 
computed directly from (17) and (18) with the expression 
for the single defect (35) and (36). Since now both A and 
A appear explicitly in the formulae for R's and T's, it is 
clear that the expansion of the double defect can not be of 
type I, but it turns out to be of type II, i.e., of the form 
(24). Hence, we will now expect that besides the even 
also the odd multiples of u> will be filtered out. This is 
confirmed by our explicit computations for two identical 
defects as depicted in figure 2. 




20 40 Q/a 60 80 

Figure 2: Absolute value squared of the Fourier transform of 
the transmission probability for a double defect with Eq — 
2.0, g = 3.5, 9 — 1.2, u — 0.2 and varying y. 

Here we have only plotted a continuous spectrum for 
y = 0.5, whereas for reasons of clarity, we only drew 



the enveloping function which connects the maxima of 
the harmonics for the remaining distances. We observe 
that now not only odd multiples of the frequency emerge 
in addition as harmonics, but also that we obtain much 
higher harmonics and the cut-off is shifted further to the 
ultraviolet. Furthermore, we see a periodic pattern in the 
enveloping function, which appears to be independent of 
y. Similar patterns were observed before in the literature, 
as for instance in the context of atomic physics described 
by a Klein-Gordon formalism (see figure 2 in [14]). 




Figure 3: Absolute value squared of the Fourier transform of 
the transmission probability for a double defect with Eq = 
2.0, g — 3.5, 9 — 1.2, for even and odd multiples of u = 0.2 
and varying y. 

With varying distance y the new structures can be 
modulated, i.e. we can control the intensity of certain 
peaks and also shift the cut-off. Clearly for a concrete 
application one would like the control mechanism to be 
as simple as possible and therefore it is interesting to in- 
vestigate precisely how the emission amplitude and the 
cut-off behave as functions of y. Unfortunately, this func- 
tion turns out to be not very simple as can be seen from 
figure 3, where we present our analysis for varying y and 
particular fixed harmonics. Nonetheless, there is a uni- 
versal shape common to all harmonics of the same type. 
As expected from the analytical expressions the overall 



G 



pattern for the odd and even harmonics differs. We con- 
firm our previous observation, namely that for small val- 
ues of the distance y, which corresponds to the limit of a 
single defect, there are no odd harmonics emerging. 

Let us now come to the main point of our analysis 
and see how this structure is reflected in the harmonic 
spectrum. The result of the evaluation of (33) is depicted 
in figure 4. 




5 10 15 Q / ffl20 25 30 35 40 

Figure 4: Harmonic emission spectrum for a single defect with 
Eo = 2.0, g = 3.5, 



1.2, oj = 0.2, A = i 



We observe a very similar spectrum as we have already 
computed for the Fourier transform of the transmission 
amplitude, which is not entirely surprising with regard 
to the expression (33). The cut-off frequencies are es- 
sentially identical. From the comparison between X and 
the enveloping function for T we deduce, that the term 
involving the transmission amplitude clearly dominates 
the spectrum. 

Let us now turn to the computation of the emission 
spectrum for a double defect system coupled to a laser 
field. We depict the results of our analysis in figure 5. 




Figure 5: Harmonic emission spectrum for a double defect 
with E = 2.0, g = 3.5, 6 = 1.2, u) = 0.2, A = 6. 

Once again, we observe a qualitatively similar be- 
haviour as for 7~, in particular the occurrence of even 



and odd order harmonics. We remark that in our ap- 
proach for larger values of y the normalization of the wave 
function becomes somewhat inaccurate and therefore the 
relative height in the intensities is not very precise. In 
principle this could be compensated, but for the conclu- 
sions we are trying to reach here that is not important. 

The important general deduction from these computa- 
tion is of course that harmonics of higher order do emerge 
in the emission spectrum of a defect system. 



C. Relativistic versus non-relativistic regime 

In the previous sections we have been working in an in- 
tensity regime which is close to the damage threshold of 
a solid, according to the experimental observations in [3] . 
This allowed us to see the maximum effect with regard 
to harmonic generation which at present might be visible 
from experiments. However, it is also interesting to in- 
vestigate situations which arc not experimentally feasible 
at present and of course lower intensity regimes. 

In order to judge in which regime we are working and 
whether there are relativistic effects, let us carry out var- 
ious limits. First of all we recall a standard estimation 
according to which the relativistic kinetic energy is close 
to the classical one when one is dealing with velocities 
v 2 <C 3/4c 2 . This is the same as saying that the kinetic 
energy is much smaller than the rest mass -Ekin <C ttiqC 2 . 
Making now a rough estimation for the system under 
consideration, we assume that the total kinetic energy is 
the one obtained from the laser field, i.e. the pondero- 
motive energy U p . We also ignore for this estimation any 
sophisticated corrections, such as possible Doppler shifts 
in the frequency, etc. Then the non-relativistic regime is 
characterized by the condition U p <C 1 . 
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Figure 6: Absolute value squared of the Fourier transform 
of the transmission probability for a single defect for various 
values of U p with g — 3.5, 6 = 1.2. 

Based on our previous observation that T and X ex- 
hibit a very similar behaviour, it will be sufficient here to 
study only the T in the different regimes, which will be 
easier than an investigation of the full emission spectrum 
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(33). From our analytic expression (23), we see that for 
a type I defect the quantity 7} becomes a function of U p , 
such that the regime will be the same when we rescale si- 
multaneously Eq and uj. Accordingly we evaluate numer- 
ically the Fourier transform (19), or equivalently compare 
against our analytical expression (23), and depict our re- 
sults in figure 6. 

We observe that when passing more and more towards 
the relativistic regime the cut-off is increased. The other 
feature one recognizes is that the modulating structure in 
the enveloping function of the harmonics becomes more 
pronounced. One should also note, in regard to (10), 
that the multipole structures might become more and 
more important in the relativistic regime. 

Let us now perform a similar computation for the dou- 
ble defect. From the expressions (26) and (27) we see 
that now 7}/ is not just a simple function of U p and 
therefore even being in the same regime the behaviour 
will be different when E and u> are rescaled. We alter 
in that case the regimes by rescaling Eq and keeping the 
frequency fixed. Our results are depicted in figure 7. 
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Figure 7: Absolute value squared of the Fourier transform of 
the transmission probability for a double defect for varying 
values of Eq with g = 3.5, 9 — 1.2, u) = 0.2. 

Similar as for the single defect we see that the cut-off is 
increased and the modulating structure in the envelop- 
ing function becomes more emphasized when we move 
towards the relativistic regime. In addition we note that 
the difference between the even and odd harmonic be- 
comes larger with increasing U p . This effect is more ex- 
treme for the low order harmonics. 

As a general observation we state that there are not 
any effects which seem to be special to the relativistic 
regime, but the transition to that regime seems to be 
rather smooth. 



V. CONCLUSIONS 

We have minimally coupled a laser field to an impu- 
rity system in a quantum wire described with the Dirac 



equation. Using a free field expansion for the free Dirac 
Fermion, we computed by means of standard potential 
scattering theory the reflection and transmission ampli- 
tudes associated to such type of defect. The amplitudes 
become functions of the defect coupling constant <?, their 
separation y and the laser field parameters Eq,uj and t. 
We employed these amplitudes in order to evaluate the 
emission spectrum of a dipole moment in such a system. 
Our findings for a single defect, taken to be the energy op- 
erator coupled minimally to the laser field, are that even 
multiples of the driving frequency uj are emitted. Inves- 
tigating a double defect system of two of such defects, 
we observe the emission of odd as well as even multiples 
of the original frequency. These features may already be 
observed qualitatively on the Fourier expansion of the 
transmission amplitude, even analytically. When carry- 
ing out the non-relativistic and extreme relativistic limit 
we do not observe any special effect, the transition seems 
to be rather smooth. 

There are various questions left for further investiga- 
tion. As an interesting application one may for instance 
compute the conductance in a similar fashion as in [10] 
and employ the laser to control it. 
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